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In a paper on the significance of tlie vector potential squared, Gubarev, Stodolsky and 
Zakliarov^ quote tfie following identity which holds for vector fields which vanish sufficiently 
rapidly at infinity that surface terms make no contribution: 



3 / ,3 , V ■ A(x)V' ■ A(x') + V X A(x) ■ V X A(x') 



They offer no proof of the identity other than to remark in a footnote that the relationship 
can be established by transforming to position space the simple vector identity (k x a)^ = 
fc^a^— (k-a)^ satisfied by a(k), the momentum-space counter part (i.e., the Fourier transform) 
of A(x). 

Recently, Stewart^ has provided a position-space derivation of Eq. ([T]) based on a 
Helmholtz separation of A(x) into its longitudinal and transverse contributions. In a com- 
ment on Stewart's paper, Duranct^ offers another position-space derivation, one that employs 
only basic operations of vector analysis and the solution to Poisson's equation for a point 
charge. 

In this note, we provide the derivation that appears to have been alluded to by Gubarev, 
Stodolsky and Zakharov. It makes use of the Fourier transform along with basic techniques 
of vector analysis and some elementary integration, but does not require the solution of any 
differential equations. As noted by Durand,^ this identity can be generalized to the case 
where the integrand of the volume integral involves the scalar product of two different vector 
fields, both of which diminish sufficiently rapidly at infinity. Since it is no more difficult to 
do so, we shall treat this more general case.- 

So, let A(x) and B(x) be two complex vector fields that are well behaved (i.e., that 
satisfy the Dirichlet conditions) and vanish sufficiently rapidly at infinity that there is no 
contribution from any surface terms. Then A(x) may be represented by the Fourier integral, 

A(x) = j d=^A; a(k) e^'^-'' (2) 

with its Fourier transform given by 

a(k) = ^/rf^xA(x)e-''^-, (3) 

and with a similar pair of relations holding between B(x) and its Fourier transform b(k). 
With these representations, the integral of interest becomes,- 

j d^x A(x) ■ B*(x) = J d^kj d^k' 8i(k) ■ b*(k') J d^xe'^^-^'^-"" (4) 

= (27r)3 / rf3^a(k) ■b*(k), (5) 



where, in obtaining the last hne, we have used the standard Fourier integral representation 
for the Dirac delta function, 
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3 , £xe'^''', 



{2n) 

along with its basic properties. 

Employing the familiar identity k x (k x a) = k(k ■ a) — (k ■ k)a, we may separate a vector 
field into its transverse and longitudinal parts, 

k(k - a) k X (k X a) 

and, with the aid of this decomposition, rewrite Eq. ([5]) as 

k(k - a) k X (k X a) 



(6) 



(i^a;A(x) ■B*(x) = {2nY d^k 



{2n?l d^k 



k^ 



k^ 



k(k ■ b* 



k X (k X b* 



[k-a)(k-b*) _^ (k X a) ■ (k X b* 



(7) 
(8) 



k^ 

where some basic vector identities have been employed in obtaining Eq. ([8]). Let us now 
transform back to coordinate space. Starting with the first term on the right hand side of 
Eq. ([H]), we have 

d^k 



fk 
1^ 



(k- a)(k- b*) 



(2^)6 y P 

d^k 



d'^xk ■ A(x)e" 



ikx / J3^f 



rfVk- B*(x')e 



(9) 



(2^)6 y p 



^ / d^xA(x) • Ve"^^'' / d^x'B*{x') ■ V'e^'^-^' (10) 
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i2nY 
1 

(2^ 



^ J d^xV ■ A(x) J rfVV ■ B*(x')e'''-("'-") (11) 

(12) 



,3 , ,3 ,V- A(x)V'-B*(x' 

d-^x / d-^x 



47r X — x' 



where, in obtaining Eq. ( |TT1) . we have integrated by parts, applied the divergence theorem, 
and assumed that the resulting surface term makes no contribution. To obtain Eq. f|T2|) . we 
have performed the integration over momentum space using spherical coordinates as follows: 



k"^ Jo J-1 ^ R Jo k R Jo X AtiR 

The second term on the right hand side of Eq. ([8]) is transformed to coordinate space in 
a similar manner. 

d^k 



d^k 1 
— (kxa).(kxb*) = ^ 



A;2 



d xk X A(x)e' 



-ikx 



rfVk X B*(x')e^''''(13) 



_i_ d=^xVxA(x)e-*'^"-y' (iVV'xB*(x')e^'^-"'(15) 

k/.3,.VxA(x).V-xB-(xO_ 

7 47r I X - x' I ^ ^ 

Combining Eqs. ([H]), (fT^ and flTBl) . we obtain our final result 

^3 A/ ^ ^ f ^3 /• ,3 ,[V-A(x)V'-B*(x') + Vx A(x)-V'xB*(x')] 

d-^xA(x) ■ B*(x) = d'^x d^x — (17) 

J J 47r I X — x' I 

recovering the Gubarev, Stodolsky and Zakharov result for the special case where the two 
fields are identical (and real). 

The same Fourier transform technique employed in this paper can be used to prove 
Helmholtz's theorem by which we mean the identity, 

demonstrating that a vector field V(x) is completely determined, apart from an additive 
constant, once its divergence and curl are specified. Indeed, the proof of either identity by 
the above method provides an excellent illustration or student exercise in the context of an 
undergraduate course in mathematical physics, combining as it does the methods of Fourier 
analysis and vector analysis. 
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^ Since the two vector fields are quite arbitrary, the use of the complex conjugate in Eq. dU and 
thereafter is simply a matter of formal convenience. But the notation does serve as a reminder 
that the identity holds for complex fields. 
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